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Abstract 

In this Letter, we show that under the assumption of high resolution, 
the quantization errors of fGn and fBm signals with uniform quantizer 
can be treated as uncorrelated white noises. 



1 Introduction 

Fractional Gaussian noise (fGn) and fractional Brownian motion (fBm) provide 
convenient ways to describe stochastic processes with long-range dependencies. 
Thus, they have received continuing interests in various fields and have many 
applications, e.g. modeling the communication networks flow and economic 
times series pQ-[3]. 

Of particular interest is the estimation of the Hurst exponent H of a fGn 
or fBm process. In practice, such estimations are usually done on the sampled 
and quantized time series. For example, texture images are often viewed as 2D 
fBm signals uniformly quantized to the — 255 scale [3] , [5] . As will be shown 
later, the quantization error might significantly affect the estimation result. 

To the best of our knowledge, no reports discussed the effect of quantization 
errors of fGn and fBm processes. In this Letter, we will show that under the 
assumption of high resolution, the quantization error can be viewed as a white 
noise added to the sampled fGn or fBm signal. 



2 The Discrete-Time fGn and fBm Signals 

There exist different kinds of discrete-time approximation for the continuous- 
time fGn and fBm processes, e.g. [B]-[5]. In this Letter, we will use the two-step 
discrete-time approximation signals denned in [S]: 
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1) First, the standard discrete-time fGn process W (ft) with Hurst exponent 
H G (0, 1) is defined as a weighted sequence of a standard Gaussian white noise 
W(n) 

W H (n) = | £*=0 ^n-i'W(*) for H + 3 (1) 
V ; [ Win) for # = § W 

where W(n) ~ JV(0, 1) and = rjfeg^fej , ft € NU{0}. 

2) Second, the discrete-time fBm process B (n) is represented as the running 
sum of W H (ft) 

B H {n) = W H (i) = £ E h "-7W( 3 ) (2) 

z=0 i=0 j'=0 

As pointed out in [ID ] . [IT ] . Eq.lJT]!-© are indeed ARFIMA processes, which 
can be rewritten in terms of the lag operator L: 

W H (n)-{ ^-L)^ H W{n) fovH^i 

W ( ' ~ \ W(n) for H = \ (6) 

and 

(1 - L)B H {n) = B H (n) - B H (ft - 1) = W H (ft) (4) 

where (1 — L) d = X^fe=o r(-d)r{k+i) L k ■ The truncated formulas in ([JJ)- (J4j) are 
equivalent to infinite formulas with W(i) = for i < 0; they are used here since 
we usually consider the case with W (i) = for i < [I]-[3j. 

For clarity, we focus on the above standard discrete-time fBm process but 
the conclusions can be easily extended to general cases. 

3 The High-resolution Quantization Errors of 
the fGn and fBm Signals 

The use of high resolution theory for error process analysis can date back to 
late 1940s [12], [13]. In [13], Bennett demonstrated that under the assump- 
tion of high resolution and smooth density of the sampled random process, the 
quantization error behaves like an additive white noise. In other words, the 
quantization error has small correlation with the signal and an approximately 
white spectrum; see also the good surveys in [2]-[TB]. In the sequel, we will 
show that this conclusion also holds for fGn and fBm signals. Our proof mainly 
uses the results in [17] . 

Suppose the original discrete-time fGn or fBm sequence S H (n) is bounded 
within [—6,6] in a finite time horizon [Q,i] and an M- level uniform quantizer 
in [—6, 6] is applied. We also assume the sample rate and the resolution of the 
quantizer are high enough. 
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As shown in [T7], by defining A = r , the normalized quantization noise 
e(n) of S H (n) can be represented as the normalized quantization noise of the 
sigma-delta modulator for S H (n): 

1 " / (M-l)gg(n) IV 
e(«)=2-<L( A + 2^ (5) 

i=0 v 7 

where (z) = a; mod 1 is the fractional part of x. 

In the following subsections, we will discuss the quantization errors of the 
fGn and fBm signals, respectively. 

3.1 fGn Signals 

For the sigma-delta modulator, we have the following useful lemmas. 
Lemma 1 \ i7| / Define an casual stable MA process x(n) 

71 

x(n) — ij){L)z{n) = ipiz(n — i) (6) 

z=0 

where z(n) is an i.i.d process having a certain distribution with smooth den- 
sity. If the regression coefficients ipi of this MA process satisfy that there exist 
rj > and infinitely many values of r such that \tpo + ... + ip r \ > i], then the 
distribution of the normalized quantization error e(n) under modulo sigma-delta 
modulation converges to the uniform distribution in s] under the assump- 
tion of high resolution. 

Lemma 2 \18f The following series iFo(a, z) is a special hypergeometric series 

ET(a + i) 

where a, z e C. 

If 1 ^ Re{a) > 0, the series converges throughout the entire unit circle 
\z\ = 1 except for the point z — 1. If Re(a) < 0, the series converges (absolutely) 
throughout the entire unit circle \z\ = 1. Particularly, the special hypergeometric 
series iFo(a, 1) converges to 0, when Re{a) < 0. 

Now we can prove the first main result of this Letter using Lemma 1 and 
Lemma 2. 



1 Lemma 1 is actually a slightly modified version of Property 3 in |17l , but it is not difficult 
to see that the proof in 1171 can be applied to Lemma 1 with little modification. 
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Theorem 1 The quantization error of fGn process with uniform quantizer is 
asymptotically uniformly distributed in [— o,h] under the assumption of high 
resolution. 



Proof 1 We will discuss the following three cases of < H < H = ^, and 
| < H < 1, respectively. From Eq.^, we know that the coefficients of an fGn 

signal are ipo,-. 



H- 



i) For H e (0, \), T(H - \) < 0. Because > H 



> 



Eoo 
i=0 



r(i+H-%) 
r(H-i)r(i+i) 



by Lemma 2. Thus we cannot directly apply Lemma 1 



here. However, the convergence speed of this series satisfies E"=o r^H-^)T(i+i) — 
for H G (0, ^) [19]. Based on these observations, we will derive the limit 
distribution of the quantization error through the limit of its characteristic func- 
tion. 

As proven in we can rewrite Eq. 0) as 

e{n) = l-\{e(n)) (8) 



where 6{n) = £™ =0 (2^1 

The corresponding characteristic function can be written as 



|*<9(n)>(27rl)| 
|£{exp [2tt± (W H {n) + ... + W a (0))]}\ 

(// -z) + ... + VK(0))]}| 



(9) 



The innermost sum in Eq.^j) can be grouped as 



£"=o h i 



H-l/2 



W(n - i) + ... + W(0) 



= h 



W{n) 



H- 



Hence 



lim $^ ri )(27r/) 



K 



1 )W{n-l) + 
; )W(0) 



(10) 



lim E 

n—too 



i=0 



3=0 



(11) 



Notice that the characteristic function of a standard Gaussian process is 
<J>ty(w) = exp(-i[j 2 ), lj e R. We have 



4> 



w 



3=0 



4 



< 



'„ * 1 



(12) 



TTie harmonic series X^^i 7 diverges; in other words, for any small positive 
number e > 0, we can always find a large enough integer n* such that 



or equivalently 



1=1 



(13) 



exp 



A J 2^ ,■ 



< e 



(14) 



From ill}) , hi 2(1 and {14% it follows that for any small e > 0, t/iere exists a 
large enough integer n* such that for all n > n* , 



|$<0„>(2tJ)| < 



l i=0 

< exp 



„ i 1 
AVI 



2tt- 



i n -. 
i=Q 



< e 



(15) 



which means 



lim ^<e„) (27rZ) 



1 ,Z = 
,Z 7^ 



(16) 



Therefore, the distribution of (0(n)} converges to the uniform distribution in 
[0, 1] and the limit distribution of e(n) is U[— 4, i]. 

ii) For H = j, we directly have Y^iLo^G.i = 1 7^ 0, so Lemma 1 applies. 

Hi) For H € (±, 1), X)^ = E^o r(g- + ^r(f+i) diverges by Lemma 2. 
It then follows from the definition of divergence that the premise of Lemma 1 
are satisfied. 



3.2 fBm Signals 

To analyze the quantization error of fBm processes, we need another lemma 
from [17]. 
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Lemma 3 J 17j Define an AR(1) process x(n) as 

x(n) — x(n — 1) = z(n) (17) 

If the input z(n) is a stationary independent increments, the normalized 
quantization noise e(n) of the modulo sigma-delta modulation converges to the 
uniform distribution in [—5, 5] and has a white spectrum under assumption of 
high resolution. 



Lemma 3 directly applies to the quantization error of fBm processes with 
H = \ (indeed, the Brownian motion). For H ^ 1/2 the techniques used for 
proving this lemma in [17] can also be adopted to characterize the quantization 
error. 



Theorem 2 The quantization noise with uniform quantizer of fBm process is 
asymptotically uniformly distributed and white under the assumption of high 
resolution. 



Proof 2 i) For H 



7j, we have 

B H (n) — B H (n — 1) = W H (n) 



(18) 



which follows directly from Lemma 3. 

ii) For H £ (0, i) U (5, 1), we will derive the limit distribution using the 
characteristic function. We can define 



e(n) = 1 - -(« 



(19) 



where 6{n) = Y^ =0 ( 



i=0 \ A 2 j 

From Eq.[lty and Eq.^, we know the regression coefficients of an fBm signal 
J2)=o('i~j)h'j 2 ■ Hence in the case of H G (j, 1), there exists r/ > 



B. 



are ip 

and infinitely many values of r such that 



Bfi 



> 



~ 4>B,r 
H-h 



H--. 



rh, 



■ h r 



> V 



(20) 



where the last inequality follows from part Hi) of the proof of Theorem 1. Now 
Lemma 1 applies and the distribution of e(n) converges to U[—h,h]. 

H-- H-- 

Similarly, we can prove the case of H 6 (0, |). Because h Q 2 > 0, h t 2 < 
for i > 0, we have 



rh 



H- 



H-i 



> 



(21) 
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for r £ N and thus 



lim r 



\4>B,0 + ■■■ + 1pB,r\ 



= lim r ^ c 

> lim, 

Again Lemma 1 applies. 



H-, 



,H-h 



H- 



H- 



= 



(22) 



From Eg. \19}) , the limit correlation R ((e(n), e(n + k)) can be written as 

R((e(n),e(n + k)) 
= | - E((6(n))) + E{(S(n))(5(n + k))} (23) 

where the limit mean E[x{nf) — linijv^oo X^^Li E(x(n)). 

Fork 0, ((5(n)}, (S(n + k))) converges in distribution to a random variable 
which is uniformly distributed in [0, l)x [0, 1); thus we have E((S(n)), (d(n + k)) = 

Jo Jo uv dudv = j. 

Fork = 0, we have E((S(n)) 2 ) = § and R({6(n)(5(n)) = |. 
By {HI), we have 

R ((e(n) , e(n + k)) = 

Thus, the normalized quantization error is white and is asymptotically un- 
corrected with the output of the quantized signal. Therefore, we prove the whole 
conclusion. 
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,k = 
,k^0 



(24) 



4 Some Simulation Results 

Fig.l shows a typical Power Spectral Density (PSD) for the quantization error 
of a ID quantized fBm signal, which indicates that the normalized quantization 
error is indeed white. 

Fig. 2 shows the PCA eigen-spectrum of the original fBm signal, the quan- 
tized fBm signal and the quantization error. According to the results given 
in [20]-[21], when the sampling data length if is a sufficiently large constant, 
the PCA eigenvalue spectrum of the auto-correlation of a ID fBm process with 
Hurst exponent H decays as a power-law 

\ k ~k-( 2H+1 \ k = l,...,K (25) 

It is also proven in [2D]-[3T] that the numerical eigen-spectrum of a white 
noise should be a straight line with slope «o ~ in the log-log scale (the 
slope is not strictly because of the finite sampling length effect). Moreover, 
when the ID fBm signal is corrupted with additive white noise and the SNR is 
large enough, the eigenvalue spectrum of the corrupted signal crossovers from 
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Figure 1: The PSD estimate for the quantization error of a ID quantized ffim 
signal, where H = 0.2, the quantization scale is A = 1. The PSD is estimated 
via periodogram method. 

o?i = — (2i? + 1) to arj. Comparing Fig. 2 to the simulation results provided in 
[20]- [21], we can see that under high resolution, the quantization error behaves 
exactly like a certain additive white noise. 
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Figure 2: PCA eigen-spectrum of the auto-correlation of a standard fBm signal, 
its corresponding quantized signal, and the quantization error, where H = 0.8, 
and the quantization scale is A = 1. 
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